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Abstract. We study integrals over Hermitian supermatrices of arbitrary size p + q, that 
are parametrized by an external field X and a source Y , of respective size m + n and p + q. 
We show that these integrals exhibit a simple topological expansion in powers of a formal 
parameter h, which can be identified with l/(p — q). The loop equation and the associated 
spectral curve are also obtained. The solutions to the loop equation are given in terms of 
the symplectic invariants introduced in [10] . The symmetry property of the latter objects 
allows us to prove a duality that relates supermatrix models in which the role of X and Y 
are interchanged. 
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1. Introduction 

In this article, we study supermatrix integrals that include an external field. We first 
prove that they can be expanded topologically. We then write their loop equations, and 
show that they are the same as for usual matrix integrals, and therefore, they have the 
same solution. In other words, super-matrix integrals' topological expansion is given by the 
symplectic invariants of |10j . As a consequence, we prove a duality which generalizes that of 
[HE]. Namely, we define (notations are explained below) 



(1.1) 

and we show that, if X and Y diagonal with complex entries, 

Z(m\n),(p\q)(X,Y) = Z(p|q),(m]n) (Y, X) (1.2) 

for all m,n,p,q > 0. Note that the matrix integrals Zr m i n \(p\ q \(X,Y) and Zr p \ q \r m \ n \(Y, X) 
share the same formal parameter h. The duality exchanges the size of the matrix with the 
number of sources, and it exchanges the external field Y of size p + q, with the sources X of 
size m + n. 



str Y 
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istrlMY n^isdet(x,-/M) 



Let us be more explicit. The ensemble H (p\q) of hermitian supermatrices of size p + q, is 
the set of matrices of the form: 

m ={ A cd) <"> 

where A and D are hermitian matrices of respective size p x p and q x q, and B and C = 
are fermionic matrices (entries are Grassmann anti-commuting variables) of respective size 
p x q and q x p. A Hermitian supermatrix can be diagonalized by an element of U(p\q), 
the supergroup of unitary transformations. For a short review of the theory of Grassmann 
algebras and supermatrices, see Appendices A and B, which are based on references [311131 [8], 

Consider P, a complex-valued function depending upon an Hermitian supermatrix M. Its 
expectation value with Gaussian measure in H{p\q) is defined by 

(P(M)) MeGU(plq) = [ dMe-^ IM ?P{IM)/ [ dMe-^< IM ? (1.4) 

JH(p\q) J 

where str is the supertrace and / is a supermatrix which ensures the convergence of the 
integral. We choose 

p q 

I = diag(C^,C"~i) (1-5) 
so that / is an element of U(p\q), the supergroup of unitary transformations. 

Our aim is to study the partition function of Gaussian supermatrix models containing 



sources xj as well as an external field Y: 



Z (mln)Mq) (X,Y)= e-^^y 2 / e ^ tr ^ J] sdetfc - Mf^M . (1.6) 

\ »=1 / M€GU(p\q) 
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In the last equation, X stands for a Hermitian supermatrix of size m + n with eigenvalues Xj 
(possibly not all distinct), so that Xi — M is understood as the matrix whose element in the 
jth. row and kth column is equal to Xi5jk — Mjk, while 

/ • \ /.s J +1 if i e {1, .. . ,m}, 

1—1 it !6{m+l,...,m + nj. 

It is worth mentioning that models involving hermitian supermatrix have been studied 
in the past. For instance, Itzykson-Zuber and character formulas for U(n\m) have been 
obtained in pj. This had been preceded, in in the beginning of the 1990s, by a few attempts 
to generalize the well known connection between conventional matrix models and quantum 
gravity in 2D [21 [6l [18]. It was soon realized however that supermatrix models could not 
provide a discrete version of supergravity. Especially, convincing arguments were given for 
the equivalence between supermatrix and matrix models when no external fields are involved 

Here we indeed prove, in the first section, that we can formally map the partition function 
of a matrix model to that of supermatrix model. It should be understood that this bijection 
remains true as long as the models are interpreted as linear combination of expectation values 
with respect to the Gaussian measure, like in Eq. f|2.35j) . and as long as all the matrices' size 
and entries are considered as parameters. When considering the models from a more gen- 
eral perspective, based on algebraic geometry, the relation between matrix and supermatrix 
problems becomes more subtle. In particular, we show that supermatrix models possess new 
critical behaviors and enjoy more symmetry. Especially, we will prove the following duality 
property given in Eq. (II. 21) . 



2. Topological expansion 

Here we show that the expectation value of powers of superstraces have a simple interpreta- 
tion in terms of ribbon graphs (also called fatgraphs). This naturally leads to the conclusion 
that the partition function for a matrix model of hermitian (p+q) x (p+q) supermatrices that 
contains an external field is similar to the partition function for the usual hermitian matrix 
model with an external field. As shown below, the formal power series expansion in fi? of the 
supermatrix integrals is of topological nature. This property will allow us, in section 13.31 to 
exploit the uniqueness of the solution to the loop equations. 

Proposition 1. Let y and Y be a N x N Hermitian matrix and a (p + q) X (p + q) Hermitian 
supermatrix, respectively. Define the map <j) as0 

cf)(tvy n ) = strY n Vn > 0. (2.1) 



Note that other routes have been followed for describing supergravity as a matrix models, such as the 
formal supereigenvalue models (see |17| for a review) and matrix model in superspace such as the Marinari- 
Parisi model (see [14] and references therein). 

n 

The function <p can be used only when the matrix y and the supermatrix Y as well as their respective sizes, 
are considered as parameters. We cannot for instance set N = 2 and y = diag(l, 1) and then apply cp. In fact, 
4> is a homomorphism that maps of the algebra of the symmetric polynomials in the eigenvalues of y to the 
subalgebra of the polynomials which symmetric in two set of eigenvalues of Y, yi, ■ ■ . , y p and y P +i, ■ ■ ■ , y P + q , 
that becomes independent of y v if y v = y P + q - In general <f> is not invertible. 
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Note in particular that the case n = corresponds to <fi(N) = p — q. Then 

nstrM n *e strMy \ =6\ ( T\ tr m n > e tv my ) |. (2.2) 

/ MeGU(p\q) \\ \ /m£GU(N)J 

Corollary 2. Let g be a natural number and p n = hstrY n , where Y is a Hermitian super- 
matrix of size p + q. Let moreover F = F/ m i n \(p\q\(X, Y) be the free energy of the supermatrix 
model whose partition function is Z = Z( m \ n \(p\ q \(X,Y); that is, F = — InZ. Then, the 
following formal power series holds 

F = Y,tf g - 2 F {g) (Pl,P2,---)- (2-3) 
g>o 

The latter results are in fact obvious reformulations of Theorem [8l which will be proved 
in the following paragraphs. 

Proposition 3. One has 

(^5 stry2 /M nn ...M inJn e strMyN 

Proof. By definition we have simply 

M M r strMY\ _ !dMe-^M)\^IMY {IM)iin ... { j M)in]n 

But a few manipulations give 

e str/Mr (JM) .. = a{j) 9 eStlIMY (2 6) 

OYj t 

so that 

(m M- ■ e stTMY \ - 

. d ^ 3 fdMe-^ stT W 2 e stlIMY 

a{jl) dY nil '" a[3n) dY jnin fdMe-b**"*)* " 

9 , . „ d I s t T MY 



°Wwr7-~°V")w— \ e / (2 - 7) 

Note that the order of the derivatives is important. We now make use of the Gaussian integral 
formula (|B.15[) and obtain: 

M ixh ■ ■ ■ M injn e" y ) = aUl)^~ • • • ^-^-e^ (2.8) 



Finally, we note that 



e^ Stry2 ^"^ Stryi = m- + ha ®Y ji (2.9) 
and the proposition follows. □ 



d h^ rV 2 d 
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We now wish to evaluate the expectation value of products of matrix elements by means of 
a simple generalization of ribbon graphs. We define a ribbon graph of order n as a sequence 
of n vertices graphically represented as half-edges | J. and ordered from left to right on an 
horizontal axis. Each half-edge j [ is labeled by a pair ij of positive integers. Moreover, a 
vertex j j. can be connected to at most one other vertex as long as the orientation of the arrows 
is respected; the connection of two half-edges produces an edge. FigQ] gives an example of a 
ribbon graph of order 6 and labeled by hji, ■ ■ ■ , ieje- 



FIGURE 1. A ribbon graph with 2 edges and 2 half-edges 




nji i2j2 13J3 UJ4: isJ5 ieje 



Any collection of independent graphs, {G±, G2, ■ ■ •}, will be written as a sum G1 + G2 + 

The concatenation (or union) of two independent graphs, G\ and G2, gives another graph 
G = G1G2. For instance, if 



Gi 



.11 



and G2 



12 h n h 



then GiG 2 



H Jl 12 J2 ^3 J3 



In order to take into account the presence of Grassmann odd variables, we need to equip 
the vertices with the following Z2— grading: The half-edge labeled by ij has degree e« + €j, 
where ej is defined by a(i) = (— l) £i , that is 



if ie{l,...,p}, 

1 if i £ {p + 1, . . . ,p 



q}- 



(2.10) 



The total degree of a graph or a subgraph G is equal to the sum of the degrees of all the 
half-edges contained in G; that is, if G can be separated into two independent subgraphs as 
G = G1G2, then degG = degGi + degG2 and the connection between half-edges doesn't 
affect their degree. 

FIGURE 2. Weight function W on the basic constituents (or propagators) of ribbon graphs. 



w 




d 



n ji 



nji i2j2 



We now introduce a weight function W on graded ribbon graphs whose values belong to a 
Grassmann algebra over C. We first set W(0) = 0, where stands for the empty graph. Then 
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we fix the weight of a half-edge or an edge as in Figj2j it is in correspondence with the matrix 
operations given in Eq. (|2.4p . Finally, we impose the distributivity and commutation rules of 
Figl3] and Figj4j The recursive application of these rules allows us to reduce the evaluation 
of a whole diagram's weight to a product of weights of edges and half-edges. 

The definition of the weight function is such that, when applying If on a ribbon graph G 
of order n, one can take any permutation ir of the half-edges and then multiply the weight 
of the new graph by the appropriate signum. In symbols, let G be a graph labeled by 
(iljl, • • • , i n jn), 7T a permutation of (1, . . . , n), 

K(hjl, ■ ■ ■ ,i n jn) = (i>ir{ii)j%(ji), • • ■ ) *ir(i„)J7r(i n ))) 

and let ttG denote the graph obtained by permuting the half-edges while keeping the links 
between the edges, then 

W{G) = sng(7r)iy(7rG), (2.11) 
where sng(7r) can be evaluated thanks to 

^lj'l ' ' ' ^Wn = Sn S( 7r )^ 7r (l)i 7 r(l) ' ' ' ^V(n)i7r(n) ■ (2-12) 

Consider for instance the graph given in Fig(TJ Then its weight is equal to 

In the case where p = 2, q = 2 and 

(k,ji, h,32, «3) J3, U, J4, k,35, k,h) = (1, 3, 2, 4, 3, 1,4, 4, 1,3, 2, 3), 
the weight simplifies to — ft I44Y23. 



FIGURE 3. Weight function W on independent subgraphs G\ and G*2- The signum 
depends on the degree of the components: e = deg G\ -deg G2 and e = (ei+ej) -deg G\. 
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Proposition 4. Let EG denote the sum of all possible ribbon graphs of degree n labeled by 
iiji, . . . , i n j n . In order words, let EG be equal to 

l\3l «2J2 ln3n ^'^J 
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FIGURE 4. Weight W associated to a crossing of ribbons. The signum depends on 
the degree of each ribbon, that is, e = (e^ + ej 1 )(ej 2 + e j2 ). 
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plus all distinct graphs obtained from the latter by connecting at least 2 half-edges and at most 
n or n — 1 of them, depending on wether n is even or odd, respectively. Then 



strA/y\ =W(m) 
I MEGU(p\q) 



Proof. Let G be an arbitrary ribbon graph and let 1 = id be the left identity operator on 
graphs, which means 1 o G=G. We introduce Aij, a noncommutative operator acting on G 
by creating a half-edge, labeled by ij, to the left of G. Obviously, the graph O n given in 
Eq. (]2.13p can be written as Ai 1 j 1 ...Ai n j n $. Let also B^^e be a noncommutative operator 
acting on G by connecting the half-edge labeled by ij to that labeled by kl, the first half- 
edge being located to the left of the second. We consider Bij as a nilpotent operator, which 
means BfjG = 0, and set B^j^G = if G doesn't contain a half-edge labeled by ij (or 
kl) or if all half-edges ij (or kl) in G are already connected. Thus, every ribbon graph 
of degree n can be uniquely written as a polynomial in operators Bij acting from the left 
on O n = Ai 1 j 1 . . -Ai n jjb. In particular, the sum EG of all graphs obtained from O n by 
connecting the half-edges in all possible ways, is given by 

EG = (1 + Bi 1 j lj i 2 j 2 )(l + Bi 1 j 1 i 3 j 3 ) ■ • • (1 + Bi 1 j 1: i n j n ) 

o (1 + B i2 j 2ti3 j 3 ) • • • (1 + B i2 j Stinjn ) • • • (1 + B in _ ljn _ 1)in j n )O n . (2.14) 

^From the properties of the connection operators B, we can simplify the latter expression as 
follows: 

^G = (1 + Bi 1 j 1} i 2 j 2 + Bi 1 j 1: i 3 j 3 . . . Bi 1 j lj i n j n ) 

o (1 + B i2j2>i3j3 + . . . + B i2hjinjn ) ■ ■ ■ (1 + B in _ ljn _ ljinjn )O n . (2.15) 

Now we consider the action weight function on the graph EG. We first note that 
O n ) = W(A iljl ...A ik _ ljk _ 1 )-W(£B 

t>k £>k 
The latter formula and Eq. (|2.15p then imply 

W(EG) = + Bi 1 j 1) i 2 j 2 + Bi l j 1 ^ 3 j 3 . . . Bi l j l) i n j n )Ai 1 j 1 

o (1 + Bi 2 j 2t i 3 j 3 + . . . + Bi 2 j 3t i n j n )Ai 2 j 2 • • • (1 + Bi n _ 1 j n _ 1 ^ n j n )Ai n j n $ s j , (2.16) 
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Using the fact that V ( ./, /i„ ; , .,. ,. ,1„ M C/ = 52i^ k B ikjk , ieje A ikjk G when G doesn't contain 
half-edges labeled by iijt such that t < k, we get 

W(EG) = W(D ilh ■ ■ ■ D inj J), D ikjk = A lkJk + B ikjk , kji A ikjk , (2.17) 

which corresponds the successive concatenation of a half-edge followed by its connection (or 
not) to all possible half-edges contained in the graph on its right. By making use of the 
commutation rules of Fig|3] and 0J one easily verifies that 

W(EG) =W(--- D ikjk D kh ■■■%) = (-l)K+%)K+%) W{- ■ ■ D kk D ikjk ■■■$) 

and 

d_ 



W(D ikjk D ik+ljk+1 ■ ■ ■ D inj J) = ( hY lkjk + aih)^— ) o W(D ik+l3k+1 • • • D inj J). 



Thus by induction, 

W(EG) = (hY ilh + a (ji)g^£j o ■ ■ ■ o (hY injn + a {j n )^—^ o 1, 

and the proposition follows from the comparison of the latter equation and Lemma El □ 

An immediate consequence of the latter result is a graded version of the Wick formula, 
which, in the context of matrix models, allows to express the correlation of 2k matrix elements 
in terms of a sum of monomials involving k correlations of 2 elements. 

In order to present the formula in a compact form, we need to introduce some more 
notations. Let a = (a±, . . . ,a k ) and b = (bi, . . . ,b k ) be two disjoint increasing sequences 
of integers in {1, . . . , 2k}. By V a , h = {(i ai j ai ,ibijbi), {ia k ja k ,ib k 3b k )} we denote a set of 
pairings of the indices iiji, ■ ■ ■ , i2kJ2k- Let 7r a b 6 S 2 A; be the permutation such that 

(ai,6i, . . .,a k ,b k ) = 7r aj6 (l,2, ...,2k- 1,2k). 

Define the signum of 7r a & by 



Corollary 5. Let (^Mi 1 j 1 ■ ■ ■ Mi n j n J denote the expectation value given by the left-hand- side 
Eq.^ZjQj with Y = 0. Then, 

( M i k jk M itj nl ) = hcj (Jk ) 5i k ,j e S jj , k . 
Moreover, if k is a nonnegative integer, 

Miiji ■ ■ ■ Mi 2k+1 j 2k+1 

[M ilh ■ ■ ■ M injn ) =Y,sM*a,b)(M iaijai M ibijbi y ■ ■ (M iakja M ibkjbk ), (2.18) 
where the sum runs over all distinct sets of pairings V a ,b of i\j\, ■ ■ ■ ,i2kj2k- 

We are ready to turn our attention to the expectation value of products of supertraces. 
From Proposition U] and the definition of the supertrace, we know that 

e-t stry2 (strM"e strMy ) = £ a(h) W(ET), (2.19) 
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where ST denotes the sum of all distinct graphs that we can get by connecting by pairs the 
half-edges of the following graph: 



4 k..l I ... 1 l-A.h 



A. ^}n-l%n ■ (2.20) 

Note that doted arrows ii* ii. have been used in order to ease the identification of 
components that share a same index. We set W (jx*" ij) = 1- A ribbon graph or a subgraph 
T is said to be of trace- type if it is labeled like the graph in Eq. (|2.20p . which means that its 
indices follow the pattern 



^k+i^k- 



The degree of a graph of trace-type depends only on the degree of the first and last labels, 
for the intermediate indices are always repeated. In order words, if T is of trace-type, 



deg 



T 



ik H 



(2.21) 



Notice that we adopt the following convention: if Tj is an empty subgraph of a trace-type 
graph, then 



ik ik 



(2.22) 



ik ik 

so that the weight of Tj is equal to one. 

In order to calculate the contribution of graphs of trace type, we need to establish a few 
rules that allow to decompose a graph in its independent parts. From the definition of 
the weight W and Eq. (|2.21|) . we see that if T\ and T 2 stand for 2 subgraphs of trace type 
(independent or not), then 



Ti 



T 2 



i. >• 

ik k it i 



(_l)0fc+ e <!)( e <! + e m)^ 



T 2 
e.. 1 . 




T x 


A, im^k / 



Moreover, 



W 




(2.23) 



(2.24) 



«fc*fe+i*fc+i 



This in turn implies that the previous subgraph has a nonzero contribution only if both T2 
and the edge linking ikik+i and has a degree equal to zero. In order words, these 

subgraphs contribute only if they are bosonic, so that we can commute all components of 
that type when evaluating the weight. By making use of equations (|2.2ip .([5j). (|2.24p and the 
fact that (-l)( e fc+^)( e fc+ e «) = (_i)(^+^) 5 one can establish the following. 

Lemma 6. Let Ti,... ,Ts be trace-type ribbon graphs that may be dependent or not. Then 
the equations given in Figures [5l OJ and hold. 

Let us consider a few examples in relation to the expectation value of strM n . First of 
all, one easily verifies that the contribution of the trace- type graph of Eq. (|2,20p is equal 
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FIGURE 5 . Commutation of trace type ribbon graphs 



a(i x )W 



To 



<j(i e )W 



^l n i e ^l 



To 



FIGURE 6 . Effect of an edge on trace-type ribbon graphs 

= ha(i k+1 )W 



2a 



Ti 



*i l k+i h n 



FIGURE 7. Effect of crossing edges on trace- type ribbon graphs 



E w 



%p , ip~\- 1 
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Ti 



T 4 



/ 



h 2 w 



T 4 



T 9 



T* 



%k U+i U i>x 



to HYi^hYiziz ■ ■ ■ KYi^, so summing over all the indices as in Eq. (|2,19p leads to str (hY) n . 
Consider next the contribution of a graph with one edge: 



«3 n 



From FigEJ we conclude if T^ is equal to a sequence of non-connected half-edges of trace-type, 
then 



E w 



r 3 



= frstr^r)^- 1 w 



Note that in the case where I = k + 1, which corresponds to graph with a closed loop labeled 
by ifc+i, the previous formula remains valid if we interpret sti(hYy~ k ~~ 1 as str 1 = p — q. We 
then conclude that the total contribution to the expectation value of strM n of all labeled 
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graphs Tki as in Eq.(|2]), with one edge linking the half-edges ikik+i and iiu+i, is 

<j{h)W{T M ) =nstr(^y) n+fe ^- 1 str(ny/- fc - 1 . 



l\,...,ln 



The ribbon graph ^5^7,36 given in Fig|8] is more complicated since the calculation of its 
contribution to the expectation of str M 8 requires the use of the three rules given in Figures 
El El and El One gets 

<r(h)W (T 15)27j36 ) = h 3 (p - q) sti(hY) 2 . 

i<ii,...,is<p+q 

The different factors in the last equation can be understood as follows: the 3 edges create a 
h?; p — q comes from the loop starting at i^; str(TiY) 2 is obtained when considering the two 
non-connected half-edges. 

FIGURE 8. Possible ribbon graph associated to (strM 8 ) 



r 



ML 




11 



13 



n 



l 5 lQ 17 l 8 



The previous discussion suggests that the weight associated to the sum over the indices 
of a trace-type ribbon graph depends only h and supertraces of the form str(hY) m , which 
includes the contribution str (hY)° = str 1 = p — q. It is indeed the case since Lema [6] implies 
that when reordering the subgraphs correctly, the only signs that appear in the total weight 
are factors of the form cr(ifc), each of them associated to an independent trace-type subgraph. 
This means that each independent subgraph produces a factor p — q if it is empty while it 
produces a factor stv(hY) n if it is made of n non-connected half-edges. Moreover, each edge 
produce a factor h. 

We can go further by defining a face of a trace-type ribbon graph as the region contained 
inside a loop, which is a closed path that follows the arrows (including the dotted arrow). 
Note that a face can contain an independent subgraph that contains non-connected half- 
edges. For instance, the graph given in Fig|8]has two faces: the first is delimited by the path 
i\ — > iq — > ^4 — > is — * i<i — > is — > ii; the second is found by following the path i% — > 17 — > 13. 
By making use of Lemma El we conclude that each face, when summing over all the indices, 
carries a weight of str(7iY) m , where m > corresponds to the number of non-connected half- 
edges. The weight of the whole graph is given by the product of the weights associated to 
the edges and the faces. Thus, if the labeled trace-type ribbon-graph T that contains n — m 
edges and F faces fi, each of them including rrii > non-connected half-edges, 

F 

a(i 1 )W(T) = h n ]JsiT(hY) mi (2.25) 
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Now define m = 

the total number of non-connected half-edges in T. 

Obviously, to any trace-type ribbon graph T with n — m edges, m half-edges, F faces 
and labeled by the integers i%,...,i n , we can associate a connected ribbon graph T, called a 
fatgraph star, with F faces, V = 1 + m vertices, and E = n edges labeled from 1 to n. Note 
that in the star T, all the points carrying the indices of T are considered as forming a unique 
vertex while the m non-connected half-edges of T are interpreted as being connected to m 
distinct vertices, so that the total number of vertices in T is 1 + m. For example, the fatgraph 
star associated to the trace-type ribbon graph of FigJB] is given in Fig. [9j It is a classical 
result that any connected graph with F faces, E edges, and V vertices, can be embedded on 
a surface of genus (at least) g in a way that guaranties that the edges don't cross (except at 
the vertices) if g is given by the following expression of the Euler characteristic [15] : 

x = V + F-E = 2-2g (2.26) 

Thus, if T is a trace-type ribbon graph with m non-connected half-edges, n — m edges, and 
F faces, 

F 

TT 1 <r(i\)W(T) = h n - m - 1 - F Ylhstr(hY) mi (2.27) 

ii,...,in 1=1 

which yields, when E = n and V = 1 + to, 

F F 

fr 1 <r(ii)W(T) = h E - v ~ F J\hstT{hY) m * = /r 2+2 f \\h^{hY) m ^ (2.28) 

il,...,i n i=l i=l 

Going back to Eq. (|2,19p . we conclude that if 

Y = hY, p mz =hstiY m \ (2.29) 

then 

e-3Wi str M n e s str =y T\ Pm . . ( 2 .30) 

\^ /MeGU(p\q) ^ 11^ V ' 

where the sum runs over all connected fatgraph stars V with one n-valent vertex, m univalent 
vertices and F faces, each face containing m, edges connected to rrii distinct vertices. 

Any correlation of product of supertraces str M n can be evaluated with the help of a simple 
generalization the combinatorial method exposed above. Indeed, Proposition [J] implies that 



e^*^ / ff istrM" fc ei- 

\ - LJ - n I MeGUMq) 

k=l 

= E E <T(ii)---<r(id)W(VT ni -'-T n *) (2.31) 

l<k<d l<i^,-..,i^<p+q 

where ST ni • • • T nd stands for the sum of all possible ribbon graphs that can be obtained by 
connecting (partially or totally) the half-edges contained in the following concatenation of d 
trace- type ribbon graphs (cf. Eq l2.20p : 

4 k.l I ... 1 1, il ... 4 k.i 1 ... 1 U4.h 

a. ilil V-iV . A. 4 4. % n d -X l n d (2.32) 
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FIGURE 9. The star fatgraph associated to the ribbon graph of Fig. [H The half- 
edges are labeled from 1 to 8. The graph contains V = 3 vertices (denoted a,b,c), 
F = 2 faces (denoted Fi,^), and E = 5 edges. As shown on the right-hand side, it 
can be embedded on compact surface of genus g = (E — F — V + 2)/2 = 1. 




In the last equation, the ith independent subgraph has valency rtj. The connection process 
may produce a minimum of and a maximum of d non-connected trace-type ribbon graphs. 
The use of Lemma [6] allows to eliminate the sum over the indices iH and leads to the following 
proposition which establish that any correlation function as an expansion in terms of (p—q)~ 2 
and supertraces of the external field Y . 

Theorem 7. Let p n = hstiY n , where Y denotes a Hermitian super-matrix or size p + q. 
Suppose that T is a star fatgraph with d labeled vertices, each of them having a valency rii, 
and define E = Yli=i n i- Suppose moreover that T contains 5 non-connected subgraphs Tj. 
Define the number of faces inT as F = X^=i fi> where fi is equal to the the number of faces 
in Tj. Finally, let mi stand for the number of non-connected half-edges in the ith face ofT 
and V = d + Yli=i m i- Then 

d F 

I TT -- str M nk e s strMy \ =e^Y h E - v ~ F TT p m . (2.33) 

\JLi/i /MeGU(p\ q ) ^ K ' 

where the sum extends over all possible (connected or not) star fatgraphs V. Moreover, V + 
F — E = 2 — 2g(T), where g(T) = Y^i=i9i an( ^ where gi denotes the minimal genus of the 
compact surface on which the graph Ti can be embedded. 

We remark that if h = (p — q) , then the polynomial po = 1. In that case, the previous 
theorem directly implies that when the external field is zero, the correlation functions of a 
Gaussian supermatrix model are equivalent to that of a Gaussian model involving a N x N 
matrix if one identifies h with N. 
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Corollary 8. Let H = (p — q)~ x and let n = Ylk n k be an even integer. Then, with the 
notation explained above, 

n- S trM"*e^ strMY \ =Yh 2 ^- 2 (2.34) 

h I M€GU(p\q) ^ 

where the sum extends over all possible (connected or not) star fatgraphs T in which the n 
half-edges are all connected by pairs. 

Let us make a last remark regarding the partition of the one-supermatrix model 

This is the usual form of the generating function considered for the enumeration of maps. 
As explained in [7j, partition functions of one-matrix models, when considering the i& as 
formal parameters (thus, not as complex numbers), are equivalent to the expectation value 
of products of characteristic polynomials in Gaussian ensembles. This is easily generalized 
to the supermatrix case. 

Proposition 9. Let S = diag(sj) be an arbitrary diagonal supermatrix of size (possibly 
infinite) p' + q' . Moreover, set 

p> q< 

tk = lstYS ~ k = ^ s ^- 7 £ s r* (2.36) 

«=1 i=p'+l 

where 7 is a complex parameter. Then, as long as the t^s are kept arbitrary, one formally 
has 

, ^ / I!?-] sdet(s 7 - - MV \ 

Proof. First, suppose that the eigenvalues of M are (Ai, . . . , A p , A p+ i, . . . , A p+(? ), so the eigen- 
values of LM are (Ai, . . . , A p , i A p+ i, . . . , i A p+(? ). Second, formally expand Ylk>i s tr M h /k 
in terms of the eigenvalues and get 

p' p p' p+q p'+q' 

U s ^ q) / \ nn( s r^n n 

y^strMM =JlLJzl i=p+U=p<+i . (2.38) 

^ k J p p +q p+q p 

n ^ {p - q) v k> - 1 > n n ^-w n n 

i=p'+l i=l j=p'+l i=p+l j=l 

We can rewrite the last equation as 

p' 

v Yl sdet ( s i - im v 

£ I str M k \ = (sdet $)•**-*) (2.39) 



exp 



k>l 



Y\ sdet( Sj - LM) 1 
j= P '+i 

and the proposition follows. □ 
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3. Loop equations 

In this section, we show that supermatrix models obey the same loop equations (Schwinger- 
Dyson equations) as usual matrix models. Since the solution of loop equations of usual matrix 
models is known order by order in H, in terms of invariants of a spectral curve [10], then this 
also gives the solution for supermatrix models. The only novelty, is that supermatrix models 
can have more general spectral curves than usual matrix models. 

3.1. Schwinger-Dyson equations. Schwinger-Dyson equations are called loop equations 
in the context of matrix models. They merely proceed from the invariance of an integral 
under changes of variables. 

Let us consider the change of variables / :Mh M' . The measure transforms according 
to the following rule [3]: 

dM' = J(M f , M)dM = sdet (j^) dM (3.1) 

J is the Berezinian, which is a generalisation of the Jacobian. Note that in the last equation, 
dM' /dM is a supermatrix of size (p + q) 2 when the size of M is p + q. For instance, if 



then 

dM' 



dM 



/da 85 dp §jr_\ 

| | | | 

! I % 1 

it it It, 

\ db* db* db* db* / 



(3.3) 



Now suppose that the transformation / is infinitesimal, that is 

f(M) = M' = M + eg(M) +0(e 2 ), (3.4) 
for f«l. By virtue of sdetM = exp(strlnM), we get 

J(M',M) = l + eK(M) + C(e 2 ) (3.5) 

where 

Lemma 10. Let A, B, and C be supermatrices of the same type than M . Suppose that 

9{M) = A l -—C. (3.7) 

x — BM 

Then we have the "splitting rule" 

K(M) = str (A l —— Slstrf ^r^C). (3.8) 

v ; V x-BM J \x-BM J y ' 

Proof. Firstly, we consider g\(M) = AMB. According to Eq. (|3.5[) . 

K^M) = str ULamb \ = Y f a^)a{j)^-{AMB) ij (3.9) 

hi 13 
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that is 



Ki(M) = s ^ j a{i)a{j)A ii B jj = str^strS. (3.10) 



ho 



Secondly, we choose g k {M) = A(BM) C and successively use K\(M): 

K k (M) = str ( —A(BM) k C U^str (ifM)^^) str {[BMf^C^j . (3.11) 
^ ' i=i 

We finally set g(M) = A(x - BM)~ 1 C = J2 k >o x ~ k ~ lA ( BM ) k C, so that 

k 

K{M) = J2 -^K k (M) = E E MA(BMYB)-^ m stT((BM)^C), (3.12) 



x^ 1 1 — ' * — ' 

k>l k>0 1=0 



which is equivalent to equation we wanted to prove. □ 
Lemma 11. Using the above notation, suppose that 

Then we have the "merging rule" 

K(M) = str [A C B ) . (3.14) 

V ' \ x-BM x-BM J K ' 

Proof. Recall Eqs. (|3.4p and (|3.5|) . Direct manipulations lead to 

K X (M) = str(AB) if gi (M) = A sir (MB). (3.15) 
and, as a consequence, 

k 

K k {M) = s £ j sti{A{BM) k - i C{BMf~ 1 B) if g k (M) = A str({BM) k C). (3.16) 
t=\ 

Now set g(M) = Astr ((x - BM)' l C). Hence 

K ^ = E ^r^( M ) = EX> (^m A (BM) k -^C(BMYB) (3.17) 
fc>l fc>0 t=o v ' 



and the lemma follows. □ 

Let us find the effect of the transformation / : M i— > M' = M + eg(M) on matrix integrals. 
We define the general expectation value of an analytic function G as 

(g(M)) m = J dMe~^ stlV ^ IM ^G(IM) I J rfMe-i*^(») (3.18) 

Note that the potential V is a rational function. The average of the identity matrix is 
proportional to 

J dM'e'i' ay f' M "> = J <fA/ e -i* v 'l' M >(l + C K(M) 

-^m, J ^-V(IM) + 0^ ) ) (3.1 9 ) 
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However, the measure is invariant under reparametrisation, which means 

Mm = ( 

Hence, the following equation must be satisfied: 



G(M)) m = {G(M')) m . (3.20) 



^K(PM)-Y J (l9(I f M)) i ^strV(M))'j =0. (3.21) 
Suppose additionally that 



lj OlVlij 

l J I M 



V(M) = str M k -str MY. (3.22) 

fc>0 

From 

d d 

strM fe = ka(i)M k ~ 1 j i and str MY = a{i)Yj h (3.23) 



dMij ()M,j 



we get, for any supermatrix A 

dMu 



S A ^dM~ Str V{M) = Str AV '^- ( 3 - 24 ) 



The substitution of the last relation into Eq. (|3.2ip leads to the following Schwinger-Dyson 
equation, also known as loop equation. 

Lemma 12. Supermatrix model satisfies the following loop equations: 

K{VM)\ = - (str I g(PM)V'(M)\ . (3.25) 
/ M a \ I M 

Notice that since the infinitesimal Berezinian K satisfies the same split and merge rule 
as usual matrix models, we already have that supermatrix models satisfy the same loop 
equations as usual matrix models. 

3.2. Loop equations for the supermatrix model in an external field. Here, we are 
interested in V(M) = stvv(M) - MY. 

Moreover, in order to close the set of loop equations, it was found in [10] that one should 
consider loop equations for the following expectation values, or more precisely their connected 
parts or (joint) cumulants (G) c . Specifically, we consider 

w( Zl ,...,z k ) = {UstT-^— \ (3.26) 

\i=l * / M 

II fi(y) J^j 1 \ 

u[x,y,zu...,Zk) = ( str — M str — ) 3.27 

\ x — My — r " Zi — M 

\ i=l I M 

^ / v'(x)-v'(M) fi(y) A 1 \ C 

p(x,y;z 1 ,...,z k ) = (sir [[str ) (3.28) 

\ y i=l I M 



"^Cumulants can be defined via (Gi(M ) . . . G„(M )) = X)tt{i n} rijG^{G , j 1 G : , 2 . . .) c , where the sum is 
over all partitions n of the set {1, . . . , n} while J = {jl,j2, ■ ■ ■} is an element of n. For instance (AB) = 
{AB) C + (A) C (B) C and (ABC) = (ABC) C + (AB) C (C) C + (AC) C {B) C + (BC) C (A) C + (A) c (B) c (C) c . 
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where n{y) = /x(y; Y) is the minimal polynomial of Y, i.e. the polynomial Tlj(y — Vj), where 
the product is over the disctinct eigenvalues of Y. 

Notice that u and p are polynomials in the variable y, and, if v'(x) is a rational fraction of 
x, then p(x, y; z%, . . . , z^j is also a rational fraction of x with the same poles, and with degree 
one less than v'(x). 

In terms of those expectation values, we have the following loop equations: 

Proposition 13. Let Y and M be Hermitian super-matrices of size p + q. In Eq. (|3.18p . set 

V(M) = v{M) - MY. (3.29) 

Then, for every set of variables J = {z±, . . . , zt}, we have the loop equation: 

k 

z, y; J U {x}) + ^2 u ( x i V'i !)w{x, J \ I) + ^ -jr- 



ll(r ... f i i i v i i , x ,..,,.,./,. / i, i x 8 ™&J\{ z j})-™(Zj,J\{ z j}) 



ICJ j=l dZj - - J 

= - (v'(x) -y)u(x,y;J) -p{x,y;J) + n(y)w(x,J) (3.30) 
which is the same loop equation as the usual 1-matrix model in an external field. 

Proof. We choose 

9(M) = /t 1 M f[ str (3.31) 

x — IM y — Y z« — M 

and apply the splitting rule of Lemma [10] to the Schwinger-Dyson equation of Lemma [T2l A 
few manipulations complete the proof. □ 

Since loop equations arise from local changes of variables, it is clear that loop equations 
may have many solutions, in fact as many as possible open integration domains in which 
matrices are integrated. Here, the integration domain is specified by defining our matrix 
integral as formal series, such that each term in the series is a polynomial moment of a 
Gaussian integral, like in Eq. (|2.35|) . 

Moreover we have seen from theorem [81 that our formal supermatrix integrals have a 
topological expansion 

oo 

w(x x , ...,x n ) = J2 tf g ~ 2+n w {9) (xi, ■ ■ ■ , x n ) (3.32) 

9=0 

It was found in |10| . that there is a unique formal power series solution of loop equations 
having such a topological expansion, and that unique solution was computed in terms of the 
"symplectic invariants" of a spectral curve. We explain below how to find the spectral curve. 

3.3. Spectral curve. The spectral curve can be found from the loop equation with J = 0; 
that is, 

Hu(x, y; x) = (v'(x) — y — hw(x)) u(x, y) — p(x, y) + fi(y)w(x) (3.33) 
and if we expand it into powers of h, to leading order we have: 

= ( v '{x) -y- w^(x)) u {0) (x, y) - p (0) (x, y) + fi(y)w {0 \x) (3.34) 
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We see this functional equation is greatly simplified if the second term disappears, i.e., if we 
choose y = v'(x) — w^(x). Moreover, in order to obtain algebraic relations, we introduce 
D(x), the denominator of v'{x). 

Corollary 14. Define two scalar functions of one and two complex variables respectively: 

y(x) =v'{x) -w i0) {x) (3.35) 
and E ex t which is a polynomial of its two variables such that 

D^x)- 1 E cxt (x,y) = (v'(x)-y)»{y)-p(°\x,y). (3.36) 
Then, the following algebraic equation holds: 

E eyct (x,y(x)) = 0. (3.37) 

The algebraic plane curve defined by this equation is called the spectral curve. 

As we said, the unique formal series solution of those loop equations of the form (|3.30p . 
having a topological expansion in h 2 , was computed in [10] in terms of the "symplectic 
invariants" of the spectral curve: 

£(x,y) = E ext {x,y) = 0. (3.38) 

It fact, to any spectral curve £, one can associate an infinite sequence of numbers J-^ 9 \£), 
g = 0, 1, 2, ... . The J-^'s are computed in terms of the spectral curve, and are residues of 
birational expressions of x and y. We refer the reader to [10] for detailed computations of 
the .F^'s. The result is that: 

oo 

hiZ = Y; fr 2g - 2 F {9) {£) (3.39) 

9=0 

The functions J-^ are called the symplectic invariants of the spectral curve because, as 
exposed in the following theorem, they remain unchanged if one changes the spectral curve 
without changing the symplectic form dx A dy. 

Theorem 15. [10] For all g > 2, the free energy !F^ 9 '(£) is invariant under the following 
transformations of the algebraic equation £(x,y) = 0: 

(1) x i ^ x , y i— > y + R(x) 

(2) x^cx, y^y/c ^ 4Q , 

(3) x i — y —x, y i— > y 

(4) x i — 2/, y i— > x 

where R is a rational function and c is a complex number. J 7 ^ is invariant under the 
four transformations if R is replaced by a polynomial P. J-^ is invariant under the four 
transformations up to an additive factor of iir/12. These transformations preserve, up to a 
sign, the symplectic form dx A dy. 

We stress that previous theorem implies that for all g > 0, !F^ 9 '{£) is invariant under the 
exchange of x and y in the £(x,y) = 0. This invariance will be exploited for proving the 
duality for the Gaussian model. 
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4. Gaussian model with sources and external fields 

4.1. Spectral curve. We now focus on the study of the partition function given in Eq. (|1.6|) . 
For this, we choose 

m+n 

V(M) = v(M) - MY, v(M) = -M 2 - ft In J| (x t - M) am > n{i) (4.1) 

i=l 

where M is a Hermitian supermatrix of size p + q. We suppose that the variables ccj's are 
the eigenvalues of a Hermitian matrix of size m + n. Amongst these eigenvalues only m' + n', 
say, are distinct. This means that 

m+n m'+n' 

sdet(z - X) = "[I (z - Xi) am "® = Yl ( z ~ x iT l ( 4 - 2 ) 

i=l i=l 

where the signed multiplicities dj are such that 

vnl m'+n! 

m = ^^aj and n = — a^. (4-3) 

i=l i=m'+l 

Similarly, 

p+<? p'+g' 

sdet(z -Y) = Yl(z- Vi)*™® = Yl( z ~ Vi) h - ( 4 - 4 ) 
i=i i=i 
Thus, the minimal polynomial of Y and the denominator of v' are respectively given by 

p'+q' m'+n' 

^(y) = u(y;Y) = Yl (V - Vi) and D{x) = (i(x; X) = JJ (x - x t ) (4.5) 
i=l j=l 

The spectral curve given in Corollary Q3] , now becomes 



£(x,y) = E ext (x,y) = JJ(a; - x { ) JJ(y - y ) ix-fi^ 

i j \ i 



a, 



where it is assumed that z G {1, . . . , m! + n/} and j G {1, . . . ,p y + q'}. 

Remark 16. This is more or less the same spectral curve as in the usual hermitian matrix 
model in an external field, with one notable difference. In the usual matrix model, the 
coefficients bj are necessarily positive integers, and here, supermatrix models allow more 
general spectral curves, where frj's can also be negative integers. 



4.2. Geometry of the spectral curve. The geometry of spectral curves of the type 14.6 
has been studied many times, and is a standard exercise of Riemann geometry [1 1|, 112] . Here, 
we only briefly summarize the main points, and we follow the same lines as [10] . 

The spectral curve £(x,y) = defines a unique compact Riemann surface. Generically, 
this algebraic Riemann surface has genus at most (m' + n'){p' + q 1 ) — 1. 
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The genus can be lower than that. In particular, if our supermatrix model is defined as 
perturbation of a Gaussian integral near M = 0, (as we did in Proposition [9]), then the 
spectral curve must have genus 0, i.e. it must be a rational spectral curve. 

But in general, the genus can be anything between and (m' + n')(p' + q') — 1. 

Notice that in Eq. (|4.6p . the coefficients / ( ^-tztm ^ \ have not been determined. Those 

coefficients cannot be determined by loop equations, they are related to the integration do- 
main for the supermatrix integral. Those coefficients are in 1-1 correspondence with the 
"filling fractions": 

e * = 7T- f V dx ( 4 - 7 ) 



2in i a- 

where Ai,i = 1, . . . , genus is a set of independent non-contractible cycles on the Riemann 
surface. 

A choice of integration domain, is equivalent to a choice of those coefficients, and thus is 
equivalent to a choice of filling fractions. 

The spectral curve is then described by this Riemann surface, and by two meromorphic 
functions x(z) and y{z) defined on that surface. Those two functions can be completely 
described by their poles and by their cycle integrals (|4.7j) . 

4.3. Poles. The spectral curve can be studied by determining the singularity structure of 
the algebraic equation. More precisely, let two complex function x and y. Let z belong to a 
the Riemann surface associated to the algebraic equation if £(x(z),y(z)) = for all z on the 
surface. 

In our case, it is clear that the algebraic equation (|4.6p has only three types of singularity, 
all being poles. First, it is obvious that the polynomial £ (x, y) diverges when both x and y 
become infinite. The equation £(x,y) = gives in that limit: 

MEi«i + Ej 6 i) , . 2 s h(m-n+p- q) 2 
y ~ x V On x ) = x hO 1/r (4.8) 

x x 

Second, we see that £ also diverges when x tends to X{. We call £j the point on the surface 
such that x(£i) = Xi and, in order to comply with the algebraic equation £(x,y) = 0, we 
must have in that limit: 

2/ ^__^ + (l) (4.9) 

Third, the polynomial goes to infinity if y approaches yi, which correspond to a diverging x. 
We thus set y(rji) = yi and we have in this limit: 

x ^J^ + (l) (4.10) 

y-yi 

The above equations yield the following characterization of x and y and, as a consequence, 
of the spectral curve. 

Lemma 17. The meromorphic functions x and y, have simple poles at the points: 



z = oo 
z = £i 
z = Vj 



x = oo, y = oo 

x = Xi, y = oo (4-11) 
x = oo, y = yj. 
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The spectral curve £(m\n),(v\q),(x,Y,e) ^ s complectly characterized by the residues 

Res ydx = h(m — n + p — q) 

Z—I-OO 

Res xdy = —h(m — n + p — q) 

z—>oo 

Res ydx = —Hai 

Res xdy = hb { (4.12) 

and the filling fractions 

e * = ^— f v dx ( 4 - 13 ) 

where we recall that a choice of integration domain, is equivalent to a choice of filling fractions. 

5. Duality 

We have determined the spectral curve, which allows in principle to construct the free 
energy via the following expansion: Then, when we have determined the spectral curve, we 
have: 

oo 

in ( Z (m\n),{p\q)( X i Y i e )) = ^^h? 9 ^ 2 •F fl .(£( m | n ) I ( p | 9 ) ) (.X-,y,e))) (5.1) 

5=0 

where Fg(£(m\n),(p\q),(X,Y,e)) denotes the symplectic invariant T g {£) defined in [10] for the 
particular spectral curve £ = £( m \n),(p\q),(X,Y,e)- 

Theorem 18. The following duality holds: 

Z (m\n),(p\q) ^> e ) = Z (p\q),(m\n) ( Y , X , -e) (5.2) 

Proof. Looking at the spectral curve characterized by data given in Lemma flTl we see that 
£(m\n),(p\q),(X,Y,e) an d £(p\q),(m\n),(Y,x,-e) are j us t obtained from one another by the exchange 
of x and y. We now from Theorem 1151 that the T g are invariant under such a transformation, 
so that the result follows from Eq. (|5.ip . □ 

The change e — > — e is just the change of orientation for the integration contours used to 
define the integrals. 

In case we are studying the perturbative supermatrix integral (small deformation of the 
Gaussian integral), we need a genus zero spectral curve. Indeed, in such instance, the pa- 
rameter z belongs to the complex plane, and x(z) and y(z) are rational fractions of z whose 
poles are fixed by the analysis performed in the previous section, i.e. by equations (|4.1ip to 
P~T2"|) . We therefore find: 

Lemma 19. For the perturbative matrix integral, the rational spectral curve can be parametrized 
as follows: 

J X(Z) = Z + hY,i y'{ Vi ){z- Vi ) /ro\ 

i y{z) -z n^i x ,^. • 



where the complex numbers £i, rji are obtained by solving x(£j) = Xi and y(rji) = y 
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This means that in the rational case, there is no filling fraction. The duality exposed in 
the introduction is thus a corollary of Theorem [18] and Lema [T9l 

We end this section with a simple application of the duality: we show that the expectation 
of a product of characteristic polynomials in a Gaussian supermatrix model formally tends 
to a generalization of the Kontsevich model (matrix Airy function) when l/(p — q) — ► 0. We 
first make use of Eq. (TB.13P , 

p+q 

Z (m\n),(p\q)( X >°) = Z (p\q),(m\n){0, X), and ^2<r(i)=p-q, 

1=1 

which yields 

_ ^ str A 2 f 

Z {mln)>m (iA,Q) = - — / dNe-2K st < IN ^i^ INA S det(-INr-^ (5.4) 

Zm,n\P') JH(m\n) 



In the last equation, A = diag(ai, . . . , a m + n ) and 2 m>n (/t) stands for the normalization coef- 
ficient defined in Eq. (|B,14p , Now, let W and X be a Hermitian supermatrix of size m + n 
and a diagonal matrix X with m + n complex entries, respectively. We also set 

h=(p-q)- 1 , IN = i+h 1/3 W, and A = 2 + H 2 / 3 X. (5.5) 

Finally, we formally expand the integrand in Eq. (|5.4|) in powers of K = (p—q)^ 1 and conclude 
that, as h — > 0, 

jyg str X 

6 _* (h , Z {m\n)Ml) ( 2 1 + 1 hl/3jC > °) = A1 (^) + 0(hV% (5.6) 



where 



and 



M(X)= f dWe^ stvW3+istlWX (5.7) 

J H(m\n) 



\{rn-n){h^-n) 

Z m ,n{K) = fo(m-n) 2 /3 z m,n{K)- (5-8) 



6. Conclusion 

We have proved that gaussian supermatrix integrals with external fields and sources, re- 
spectively denoted by Y = diag(yi, . . . y p + q ) and X = diag(xi, . . . , x m+n ) for some complex 
numbers yi and Z{ , satisfy a duality formula which extends that of [U [7] . Usual hermitian 
matrices with sources at both numerators and denominators could not have this duality, 
because we see that numerators and denominators are transformed, under this duality, into 
variables of different signs, which can be only obtained with supermatrices. In some sense 
supermatrices allow eigenvalues with negative multiplicities. We conjecture that the super- 
matrix duality proved in the article extends to the case where the external sources and fields, 
X and Y, are arbitrary supermatrices. 
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Appendix A. Grassmann variables 

The Grassmann algebra of order n, denoted T = T(n; C), is the algebra over C generated 
by 1 and n quantities 0j, called Grassmann variables, which satisfy 

6i0j = -Ofr. (A.l) 

Note in particular that Of = 0. Let c be a complex number. We adopt the following convention 
for the complex conjugation: 

(c0i)* = c*9t, (0 i j y = 9*9*, 9** = 9 l . (A.2) 
T has 2 n generators: Every element x can be decomposed as 

x = x (o) + x^9 h + x {il ' i2) 9 h 9 i2 + ... + x^'-riO! ■ ■ ■ 9 n (A.3) 

il i\<t2 

with ar(*i'*2'-) g c \y e sa y that x is even (or bosonic) if it contains only monomials with an 
even number of 9f, x is odd (fermionic) if contains only monomials with an odd number of 
9j. Elements the Grassmann algebra satisfy 

x iXj = (-l)^ Xi - de ^x jXi (A.4) 

where deg x = if x is even while deg x = 1 if x is odd. Note that we also write 

degx = e(x) and (-1)*** = a (x). (A.5) 

The (left) derivative and the integration with respect to the Grassmann variables are 
respectively defined by 

-QQ^h^h " " " Ojk) = S hh (9j2 ■■■Ojk)- S i,32 (^1^3 ■ ■ ■ 9jk) + ( -1 ) fc ^,J* (Oji • • • %fc-i) ( A - 6 ) 

and 

J dOi 9j = - J 9j dOi = Si j. (A. 7) 

With the latter definition, one can obtain the following integral representation for the deter- 
minant of a N x N matrix X: 

detX= I d9U9e- 0tX0 



J d9U9e- e ' xe (A.8) 
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where (resp. 9*) is interpreted as a column (resp. row) vector with ./V fermionic components. 
The bosonic counterpart of the latter formula is the usual Gaussian integral involving complex 
variables zf, that is, if X is Hermitian, 

i i r . + . _, tx 



dzUze-* x *. (AS 



det X (2vr i)^ 

Appendix B. Supermatrices 
Let X denote an even (bosonic) square supermatrix of size (p + q). It can be written as 

where A and D are respectively p x p and q x q matrices with even Grassmann elements, 
while B and C are respectively p x q and gxp matrices with Grassmann odd elements. The 
following notation is useful: 

(-l)** x v=a(i)a(j) (B.2) 

where 

= (-!)««) = i=1 -"> (B.3) 

[-1, i = p+l,...p + g. 

The supertrace is given by 

p+g 

sir X = a(i)Xu = tvA-trD (B.4) 

i=i 

and satisfies 

str (X + Y ) = str X + str Y, str(XY) = str (YX) . (B.5) 
The superdeterminant, which exists only if A as well as D are invertible, is given by 

, v det( A- BD^C) det (A) s 

Sdet * = det( D ) - det(D- cI-B) - (R6) 

One can show that 

sdet (XY ) = sdet X sdet Y and det(expX) = exp(strX). (B.7) 
We use the following definitions for the transpose and the adjoint of a supermatrix: 

*-(* g) and Xt-(^ (B.8) 

where ^ means (A*)*. On easily shows that 

(X^ 1 " = X, (XY) f = Y f X f , sdetX* = (sdet X)*. (B.9) 

However, (XY)' 7^ Y t X t and sdetX* 7^ sdetX in general. A supermatrix X is Hermitian 
if X^ = X; it is unitary if X^ = X -1 . The set of all invertible even supermatrices of size 
p + q, whose elements belong to a Grassman algebra over C, form the general linear super 
group GL(p\q). All unitary supermatrices form the superunitary group U(p\q). 



4 Other non equivalent definitions for the transpose and the adjoint are possible. They lead, for instance, 
to a distinct unitary supergroup, namely sU(p\q). [51 113]. 
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Let £' = (zl, . . . ,z*,9*, . . . , 6*) be the adjoint of the supervector £. Then, the superdeter- 
minant of a Hermitian supermatrix X has the following integral representation: 

-i- = -V / det^e"^^. (B.10) 
sdetX (2vri)P J v ; 

The measure for Hermitian supermatrices is given by the product ■ dXij which is equal, 
up to a multiplicative constant, to the product of real independent differential elements. 
Explicitly, if X G H(p\q), 

dX= Yl dX a n dXijdXfj jj dXijdX*j. (B.ll) 

l<i<p+q l<«<i<P 1<«<P 

P+l<i<p+<? 

For instance, by using the bosonic and fermionic Gaussian integrals, 

[dxe- x2+xy = v^Fe^ 2 and f «*e ww = e™', (B.12) 

where degry = degry* = 1, one readily shows that the following formulas hold for any (not 
necessarily Hermitian) supermatrix Y of size p + q: 

J (|Ie^ str ^ 2 e str/1Y = z Piq (h)e% stvY \ (B.13) 

where 

Z p , q (H) = 2 { P+^' 2 iM vr (P 2 +9 2 )/2 n ( P -9) 2 /2 ) (B . M) 

so that 

/ p strXr\ J dle a e _ r SstrY 2 m 15^ 
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